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Die Erfindung betrifft ein Verfahren zur Bestimmung der, insbesondere zukunftigen, 
Position eines Objekts, insbesondere eines Teilchens Oder dergleichen, wobei der 
Positionsverlauf eines Referenzobjekts bestimmt wird, wobei das Referenzobjekt 
einen Prozess durchlauft und wobei die Differenzen der Positionen des Objekts und 
des Differenzobjekts gebildet werden. 

In vielen physikalischen und in verschiedenen technischen Systemen ist es notwenig 
die Position, insbesondere die zukunftige Position, eines Objekts zu kennen, um bei- 
spielsweise Schaden zu vermeiden, die dadurch entstehen, class das Objekt sich in 
einer Position befindet, in der ein solcher Schaden verursacht wird. Zur Bestimmung 
der Position sind hierbei verschiedene Verfahren bekannt, bei denen der Positions- 
verlauf eines Referenzobjekts bestimmt und die Differenzen der Positionen des Ob- 
jekts und des Referenzobjekts gebildet werden. Bei dem Referenzobjekt handelt es 
sich hierbei um ein Referenzobjekt, dessen Dynamik gut bestimmbar ist. Bei den be- 
kannten Verfahren ist insgesamt problematisch, dass diese oftmals nur unzureichen- 
de Ergebnisse liefern, da die zur Bestimmung der Position verwendeten Methoden 
nur mit einem technisch sehr hohem Aufwand genaue Ergebnisse liefern. Steht die- 
ser technische Aufwand nicht zu Verfugung, so sind die Ergebnisse oftmals nicht 
* brauchbar. 

Der vorliegenden Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren zur Be- 
stimmung der Position eines Objekts der eingangs genannten Art anzugeben, bei 
dem eine gute Bestimmung der Position des Objekts bei gleichzeitig moglichst gerin- 
gem technischen Aufwand ermoglicht ist, 

ErfindungsgemaR wird die voranstehende Aufgabe durch das Verfahren zur Bestim- 
mung der Position eines Objekts mit den Merkmalen des Patentanspruchs 1 gelost. 
Danach ist das in Rede stehende Verfahren zur Bestimmung der Position eines Ob- 
jekts derart ausgestaltet und weitergebildet, dass aus den Differenzen Fehlpositions- 
terme gebildet werden und dass die Fehlpositionsterme mit mindestens einem Kor- 
rekturfaktor gewichtet werden. 

In erfindungsgemaRer Weise ist erkannt worden, dass durch die Einfuhrung von ge- 
wichteten Fehlpositionstermen in besonders einfacher Art und Weise der technische 



Aufwand, so beispielsweise der Rechenaufwand bei Simulationen, gering gehalten 
werden kann. 

Hinsichtlich einer besonders einfachen Ausgestaltung konnten die mit dem Korrektur- 
faktor gewichteten Fehlpositionsterme aufaddiert werden. Somit ware ein Korrektur- 
. faktor gegeben, der besonders einfach ermittelbar ist und au3erdem eine gute Kor- 
rekturleistung aufweist. Der Korrekturfaktor konnte in weiter vorteilhafter Weise durch 
eineh adjungierten Prozess berechnet werden. Bei dem adjungierten Prozess konnte 
es sich hierbei um einen Prozess handeln, der in entgegengesetzter Zeitrichtung des 
Prozesses des Referenzobjekts verlauft. Im Rahmen einer besonders einfachen 
Ausgestaltung konnten zur Bestimmung der Position des Objekts die gewichteten 
und aufsummierten Fehlpositionsterme zu der Position des Referenzobjekts addiert 
werden. 

Im Rahmen einer besonders einfachen Ausgestaltung konnte es sich bei dem Objekt 
auch um eine amerikanische Option und bei dem Referenzobjekt um die aquivalente 
europaische Option handeln. Bei der Position konnte es sich dann um den fruhen 
Ausubungspreis, bei dem Fehlpositionsterm wurde es sich in diesem Fall um den 
Fehlpreisterm und bei dem Korrekturfaktor um die Konstante handeln. 

In erfindungsgemaBer Weise konnte dann fur eine vorgegebene Optimierungsme- 
thode, beispielsweise der Monte-Carlo- oder der Finite-Differenz-Methode, der fruhe 
Ausubungspreis berechnet werden, das heiGt die Preisdifferenz zwischen amerikani- 
schen und europaischen Optionen. Vorzugsweise konnte dies durch die Summierung 
eines Fehlpreisterms erreicht werden, der bereits durch die meisten Methoden be- 
rechnet wird. Dieser konnte in besonders vorteilhafter Weise durch eine Konstante 
gewichtet werden, die durch einen entsprechend ausgewahlten adjungierten Prozess 
berechnet werden konnte. Der adjungierte Prozess konnte sich in weiter vorteilhafter 
Weise in entgegengesetzter Zeitrichtung des ursprunglichen Prozesses entwickeln, 
der die Preisentwicklung der europaischen Optionen bestimmt, wobei beispielsweise 
ein Entwicklungsparameter der formale Verknupfungsoperator des europaischen 
. Preisoperators sein konnte. Durch die Schatzung der Pramie fur die fruhzeitige Aus- 
Cibung - Early Exercise Premium - konnte dann der Preis von amerikanischen Opti- 
onen einfach dadurch berechnet werden, dass er zum Preis der entsprechenden eu- 



ropaischen Option addiert werden wurde. Hierdurch konnten somit numerische Feh- 
ler in den mangelhaften Simulationsmethoden kompensiert werden und bessere Ab- 
schatzungen erreicht werden als bei den reinen numerischen Methoden, auf die das 
erfindungsgema(3e Verfahren aufgesetzt wird. Dies erlaubt beispielsweise eine 
schnelle Simulation basierend auf den groben Gittern der Finite-Differenz-Methode. 
Es ist so auch moglich praktische Grenzen der zugrundeliegenden Methode zu u- 
berwinden und macht die Bewertung von komplexeren Optionen mit einer Vielzahl 
von Basiswerten durchfuhrbar. 

Das Verfahren ist insbesondere ftir die Schatzung von amerikanischen Standardopti- 
onen mit einem einzigen Basiswert geeignet und besonders vorteilhaft einsetzbar als 
Zusatz zur Finite-Differenz-Methode und der Monte-Carlo-Methode. Die Realisierung 
erfolgt in weiter vorteilhafter Weise in Java. 

Das Verfahren konnte auch auf verschiedenste, insbesondere regelungstechnische, 
gesteuerte oder ungesteuerte Prozesse anwendbar sein. In Bezug auf die bereits 
erlautere Ausgestaltung bezuglich der Schatzung von Optionspreisen konnte der 
Preis einer europaischen Option ein ungesteuerter Prozess sein, da er unabhangig 
von auBeren Einflussen des Inhabers ist. Bei der Bestimmung des Preises einer a- 
merikanischen Option konnte es sich entsprechend urn einen gesteuerten Prozess 
handeln, da in diesem Fall eine optimale Preisbestimmung erfolgen soli, d. h. das die 
Option ausgeubt werden konnte, sobald der Preis unter den inneren Wert der Option 
fallt. So lasseri sich auch die weiteren Ausfuhrungen ohne weiteres auf gesteuerte 
und ungesteuerte Prozesse anwenden und ubertragen. 

Es gibt nun verschiedene Moglichkeiten, die Lehre der vorliegenden Erfindung in vor- 
teilhafter Weise auszugestalten und weiterzubilden. Dazu ist einerseits auf die dem 
Patentanspruch 1 nachgeordneten Patentanspriiche und andererseits auf die nach- 
folgende Erlauterung bevorzugter Ausfuhrungsbeispiele des erfindungsgemaf3en 
Verfahrens zur Bestimmung der Position eines Objekts anhand der Zeichnung zu 
verweisen. In Verbindung mit der Erlauterung der bevorzugten Ausfuhrungsbeispiele 
des erfindungsgema3en Verfahrens zur Bestimmung der Position eines Objekts an- 
hand der Zeichnung werden auch im Allgemeinen bevorzugte Ausgestaltungen und 
Weiterbildungen der Lehre erlautert. In der Zeichnung zeigt 



Fig. 1 ein Ausfuhrungsbeispiel der Abhangigkeit einer amerikanischen Option 

bei einem risikolosen Zinssatz r= 0,1, der Volatilitat a= 0,4, dem Aus- 
ubungspreis E = $ 10, wobei die Option in sechs Monaten ablauft, 

Fig. 2 die Diskretisierung der Preise der Basiswerte sowie die zugelassen 

Ubergange fur einen Zustand und 

Fig. 3 ein weiteres Ausfuhrungsbeispiel mit der Darstellung des Preises einer 

amerikanischen Put-Option auf einem groben Gitter mittels einer klei- 
nen Anzahl von Ubergangen (MC30) und durch Hinzufugen der adjun- 
gierten Prozesskorrektur in der gleichen Simulation des geschatzten 
europaischen Optionspreises (APC) mit einer Volatilitat a - 0,4, einem 
Zinssatz r= 0,1 und dem Ausubungspreis E= $10. 

Artierikanische Optionen unterscheiden sich von europaischen Optionen dadurch, 
dass sie jederzeit vor dem Ablauf der Option ausgeubt werden konnen, wahrenddes- 
sen europaische Optionen nur an ihrem Falligkeitsdatum ausgeubt werden konnen. 
Das Schatzen von amerikanischen Optionen ist daher weitaus schwieriger als das 
Schatzen von europaischen Optionen und erfordert speicher- und zeitintensive nu- 
merische Simulationen. 

Es ist bekannt, dass amerikanische and europaischen Optionen unter bestimmten 
Voraussetzuhgen Equivalent sein konnen, beispielsweise fur eine amerikanische 
Standard-Put-Option auf eine Aktie, fur die keine Dividende ausgezahlt wird. Im All- 
gemejnen muss zur Berechnung eines angemessenen Preises fur eine amerikarii- 
sche Option in Betracht eine optimale Ausubungsstrategie zu jedem Zeitschritt gezo- 
gen werden, was analytische und numerische Berechnuhgen erheblich umfangrei- 
cher macht als im Falle der europaischen Optionen. 

Die Schatzung von amerikanischen Optionen ist somit komplexer als die Schatzung 
von europaischen Optionen. In dieser Arbeit wird eine Formel entwickelt, die es er- 
laubt einen fruhen Ausubungsoptionspreis berechnen, d.h. die Preisdifferenz zwi- 
schen den beiden Optionstypen im Hinblick auf einen adjungierten Prozess, der sich 



in der sich in der ruckwartigen Zeitrichtung des ursprunglichen Prozesses entwickelt, 
der die Entwicklung des Preises der Europaischen Optionen bestimmt. Diese Be- 
rechnung kann benutzt werden, urn die Optionspreisschatzungen zu verbessern, die 
durch numerische Methoden, wie beispielsweise die der Finiten-Differenz- oder der 
Monte-Carlo-Methode, erhalten werden. 

Es wird ein Verfahren entwickelt, das als Adjungierte Prozesskorrektur Methode - 
Adjoint Process Correction Method (APC) - bezeichnet wird, welches die Berech- 
nung des Preises einer amerikanischen Option durch die Korrektur des Preises sei- 
nes europaischen Aquivalent erlaubt. Es wird seine Zweckmaftigkeit durch die An- 
wendung der APC-Methode auf zwei verschiedene numerische Methoden zur Opti- 
onspreisberechung gezeigt, namlich die der Finiten-Differenzen und der stochasti- 
schen Simulation. Die letztere basiert auf einer Umformulierung einer Entwicklungs- 
gleichung im Hinblick auf eine Mastergleichung und wird haufig in physikalischen 
Prozessen verwendet, urn eine Vielzahl von sowohl stochastischen als auch deter- 
ministischen Phanomenen zu modellieren und zu simulieren. Anders als die Finite- 
Differenz-Methode kann die stochastische Simulation fur praktisch jede bekannte 
Preisdynamik verwendet werden. 

Die APC-Methode verwendet hierbei sowohl die bekannte Ruckwartsentwicklungs- 
gleichuhg als auch einen adjungierten Prozess, der vorwarts in Zeitrichtung ablauft. 
Die Preiskorrektur konnte somit als ein Integral iiber einen Quellterm ausgedruckt 
werden, der die Preisdiskrepanz zwischen den europaischen und den amerikani- 
schen Ruckwartsentwicklungsoperatoren fur infinitesimale Zeitschritte beschreibt, die 
durch den adjungierten Prozess gewichtet sind. Oftmals konnen die notwendigen 
Informationen uber den adjungierten Prozess analytisch bestimmt werden. Wenn 
dies moglich ist, kann die APC-Methode besonders einfach und bei gleichzeitig ge- 
ringfugigen Kosten als eine Erweiterung von konventionellen Methoden benutzt wer- 
den. In diesen Fallen kann die APC-Methode Informationen genau entwickeln, die 
anderenfalls durch Diskretisierung oder statistische Fehler gestort wurden und kann 
somit bei nur geringfugig hoherem Rechenaufwand genauere Ergebnisse als die ein- 
fache Finite-Differenz-Methode oder die Monte-Carlo-Methode liefern. 



Ein weiterer Vorteil der APC-Methode ist, dass sie auf die relevanten Preiszonen fo- 
kussiert werden kann. Wenn beispielsweise der derzeitige Preis eines Optionsbasis- 
werts $ 10 ist, ist der Optionspreis fur $ 20 nicht besonders interessant - aber im All- 
gemeinen wird die Finite-Differenz-Methode diesen Wert trotzdem angeben, genauso 
wie alle anderen Werte auf einem Feingitter zwischen $ 10 bis $ 20, da das $ 20- 
Gebiet wahrend der Ruckwartsrechnung der Endbedingungen von entscheidender 
Wichtigkeit sein konnte. Der adjungierte Prozess ermoglicht es, zwischen $ 10 und $ 
20 ein vergleichsweise grobes Gitter zu verwenden und trotzdem sehr genaue Er- 
gebnisse fur den derzeitigen Preis durch die Verdichtung der Anfangsverteilung des 
adjungierten Prozesses bei $ 10 zu erhalten. Fur Optionen mit mehrfachen Basiswer- 
ten wachst die Komplexitat des Problems exponentiell mit der Anzahl der Basiswerte, 
demnach muss fur d Basiswerte ein tf-dimensionaler Raum diskretisiert werden. 
Deshalb ist die Fahigkeit der APC Methode, auch bei groben Gittern gute Ergebnisse 
zu liefern, besonders vorteilhaft bei numerisch anspruchsvollen, hoher dimensionalen 
Fallen. 

Im numerischen Teil der Arbeit werden Ergebnisse fur einen einzigen Basiswert ge- 
zeigt, der einer geometrischen Brown'schen Bewegung folgt, und wird mit dem 
Black-Scholes-Modell in der selbst-adjungierten Form gearbeitet. Hierbei ergeben 
sich allerdings keine Einschrankungen der APC-Methode an sich. Beispielsweise 
kann die stochastische Simulation auf komplizierte Prozesse generalisiert werden 
und Sprungprozesse konnen in einfacher Art und Weise integriert werden. AuBerdem 
besteht keine Notwendigkeit fur den Zeitentwicklungsoperator selbst-adjungiert zu 
sein. Im weiteren wird von einem stetigen Zeitverlauf ausgegangen und angenom- 
men, dass die Entwicklung des europaischen Optionspreisvektors P wird durch die 
Gleichung 

bestimmt wird, wobei die APC Methode in einfacher Art und Weise auf diskrete 
Zeitverlauf e ubertragen werden kann. 

Zunachst wird gezeigt, wie der Preis von einer europaischen Option in der Art eines 
adjungierten Prozesses ausgedruckt werden kann. Dieses Ergebnis wird benutzt und 



generalisiert auf den Fall einer amerikanischen Option, was zu der Korrekturformel 
fuhrt, die in der APC-Methode benutzt wird. In einer ersten numerischen Anwendung 
wird gezeigt, das die APC-Methode dazu verwendet werden kann, urn Optionspreise 
durch die Finite-Differenz-Methode mit einem sehr groben Gitter zu berechnen und 
zu interpolieren. Eine alternative stochastische Methode wird dargestellt, die zeigt, 
dass die Optionspreisberechnung durch den adjungierten Prozess eine sinnvolle 
Preisschatzung ergeben kann, wenn die Statistik der zugrundeliegenden Monte- 
Carlo-Methode zu schlechte Ergebnisse liefert, um sinnvolle Ergebnisse zu liefern. 

Zu Beginn wird der Preis einer europaischen Option in der Art eines adjungierten 
Prozesses entwickelt sowie Definitionen und invariante Eigenschaften erlautert. Die- 
ser Prozess wird zunachst definiert und eine erhaltene Quantitat gefunden. 

Bei einen adjungierten Prozess des Prozess Pder Gleichung (1) ist ein Prozess ge- 
meint, der durch die Gleichung 

definiert wird, wobei L* die formal Adjungierte des Operators mit dem Standard L 2 
Skalarprodukt, so dass 

fur alle y/, </> im Definitionsbereich von L beziehungsweise L* gilt. Es wird angenom- 
rrien, dass alle (f> in dem Definitionsbereich L* stark abfallende Funktionen sind oder 
fur einen endlichen Definitionsbereich gering im Grenzbereich verschwindend sind, 
um Probleme im Randbereich bei der partiellen Integration zu vermeiden. Insbeson- 
dere wird dies als Randbedingung fur P* vorausgesetzt. Wenn beispielsweise die 
Dynamik durch eine Advektionsdiffusionsgleichung durch L = UV + DA bestimmt 
wird, ist der adjungierte Operator L* = - UV + DA, woraus ersehen werden kann, 
dass der Prozess selbst in der positiven entwickelnden Zeitrichtung eindeutig defi- 
niert ist, wahrend der adjungierte Prozess in ruckwarts entwickelnden Zeitrichtung 
eindeutig definiert ist. In der sich vorwarts entwickelnden Zeitrichtung ist die Losung 
der adjungierten Prozessgleichung aufgrund der Anwesenheit des negativen Diffusi- 




onsoperators eine schlecht losbare Aufgabe. Entsprechend ist die zeitliche Zuruck- 
entwicklung fur den ursprunglichen Prozess schlecht losbar. 



Mit diesen Definitionen ist das Skalarprodukt von P und P*zeitlich invariant. Tatsach- 
lich gilt wegen der Definition des adjungierten Operators in Gleichung (3) fur irgend- 
welche Zeitpunkte T 1 und T 2 \ 



Die Anfangs- oder Endbedingungen fur den adjungierten Prozess konnen nun in ei- 
ner Art gewahlt werden, die am Geeignetesten fur das zu losende Problem ist. Zu 
diesem Zweck wird 



gewahlt, wobei x 0 der Assetpreis ist, an dem die Option bewertet werden soli und 8 
die tf-dimensionale Dirac-Distribution ist. In vorteilhafter Weise konnte x 0 als ein zu- 
satzlicher Parameter eingefuhrt werden und konnte durch P*(x,t; x 0 ) die Losung des 
Anfangswertproblems der Gleichungen (2) und (5) bezeichnet werden. Mit diesen 
Anfangsbedingungen ergibt sich 




0 



(4) 



(5) 




(6) 



wobei in Gleichung (4) T 1 = 0 und T 2 = Tgesetzt wird, ohne das die Allgemeinheit der 
Formel beeintrachtig wird. 



Die Schatzung der Option erfolgt nunmehr mittels eines adjungierten Prozesses. Es 
wird vorausgesetzt, dass der Preis einer europaischen Option durch folgende Glei- 
chung bestimmt wird: 

(7) 

In dem Modell von Black & Scholes und Merton ist beispielsweise L = L BS der Black- 
Scholes-Operator fur einen risikolosen Zinssatz r und einer zugrundeliegenden Vola- 
tility <r, der gegeben ist durch 

£bs'?= -o^*m :t ~r&m ( 8 ) 



Aus dem Vorzeichen des Diffusionsteils - Vz a 2 S 2 d 2 /d S 2 des Operators ist ersicht- 
lich, dass die Gleichung riickwarts in der Zeit gerechnet wird. Die Endbedingungen 
fCireine Put-Option an ihrem Falligkeitstag 7 ist 

tl =-ma*(i^ mfymwiM,. (9) 

wobei E der Basispreis der Option ist. Mit dem entsprechenden adjungierten Prozess 
P* der durch 

p*&Mfr stx -$) (10) 

gegeben ist und Gleichung (6) erhalt man 

- f dx P*\x,T\ S) max(rc - M, 0) 

J ' (11) 



wodurch der Preis einer Option zum gewartigen Zeitpunkt t = 0 und gewartigen Preis 
S in den zugrundeliegenden Bedingungen des adjungierten Prozesses ausgedruckt 



werden kann. Im Black-Scholes-Modell entspricht dies der Optionspreisberechnungs- 
formel, die den Preis als diskontierten Erwartungswert eines Profits bei einer risiko- 
neutralen MaBgabe bestimmt. 

Der adjungierte Prozess wird nunmehr auf die Schatzung von amerikanischen Optio- 
nen angewandt. Fur amerikanische Optionen beinhaltet die Preisfindung bei jedem 
Zeitschritt eine Entscheidung, ob die Option ausgeubt werden soil oder nicht. Dies 
macht die Dynamik des Prozesses komplizierter und verschiedene Herangehenswei- 
sen sind hergeleitet worden, um diese Dynamik zu modellieren, beispielsweise als 
Free-Boundary-, Linear-Complementarity- oder als Optimal-Stopping-Problem. Hier- 
bei wird angenommen, dass fur infinitesimale Zeitschritte dt die Zeitentwicklung des 
Optionspreises durch folgende Gleichung dargestellt wird: 



wobei F(x, t) ist eine Payoff-Funktion, die den intrinsischen Wert der Option wieder- 
spiegelt, der hier sowohl vom derzeitigen Status xdes Basiswerts als auch der Zeit t 
abhangt. Diese Formulierung kann beispielsweise aus dem Black-Scholes-Modell als 
eine explizierte Diskretisierung des entsprechenden linearen Komplementarprob- 
lerns, wie von Wilmott, Howison und Dewynne gezeigt, hergeleitet werden. Ahnliche 
Gleichungen konnen auch aus diskreten Zeitformulierungen, beispielsweise der Dy- 
namic-Programming-Gleichung fur binomische Baume, hergeleitet werden. 

Unter der Berucksichtigung von 



time (Fan + £$a^d%$) ^#am + + - I^m ^ 0} (13) 



(12) 



kann die Gleichung fur dt} 0 umgeschrieben werden zu 




(14) 



mit einem Quellterm 



- 12 - 



Dabei ist zu beachten, dass der Quellenterm nur in dem Gebiet der Ausubung nicht 
verschwindet. In dem Gebiet der Ausubung gilt jedenfalls P am (x, t) = F(x, t), wodurch 
folgendes erhalten wird: 

0 i (16) 

Aus den Gleichungen (15) und (16) konnte auch ohne numerische Simulationen be- 
stimmt werden, ob die amerikanische Option sich von der europaischen Option un- 
terscheidet. Dies ware insbesondere dann moglich, wenn gilt, dass (F - P am - 
LP am df) und/oder (d/dtF-LF) < 0 ist. In diesem Fall verschwindet der Quellterm J(F f 
P a m> 0. woraus folgt, dass sich die beiden Systeme im Wesentlichen identisch verhal- 
ten. Ein Ausfuhrungsbeispiel fur diesen Fall ist eine amerikanische Standard-Put- 
Option auf Aktien ohne Dividende im Black-Scholes-Modell. 

Hierbei ist besonders vorteilhaft, dass numerischer Aufwand vermieden werden 
konnte, da man nunmehr den leichter zu ermittelnden Preis der europaischen Option 
ermitteln konnte. AuBerdem konnten somit gezielt Produkte geschaffen werden, die 
sich im Wesentlichen wie eine europaische Oder amerikanische Option verhalten. 
Durch die Wahl der Payoff-Funktion F lasst sich das Verhalten der Option derart be- 
einflussen, dass J(F, P am , t) nicht negativ oder nicht positiv wird. Eine fruhzeitige 
Ausubungsmoglichkeit konnte fur beide Seiten insbesondere dann vorteilhaft sein, 
wenn dem Anbieter bekannt ist, ob eine Option ohne oder mit nur geringem Risiko- 
aufschlag als amerikanisch angeboten werden kann. 

Es wird nun der adjungierte Prozess P*wie in Gleichung (2) mit L der Gleichung (12) 
definiert, das heiBt anders als P ami wird der adjungierte Prozess wie bei einer europa- 
ischen Dynamik bestimmt. Wird die Berechnung der europaischen Optionen wieder- 
holt und wird angenommen, dass die derzeitige Zeit Null und die Ausubungszeit 7 ist, 
so gilt ohne Beschrankung der Allgemeinheit, 
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j dx (p^x, T)P'(x f T) - P am {x, 0)P>,0)) 

= £ di j dx—firrfc^P'fa *)) 

= jT# J dx(p%x,t)CP arn {x,t) - *UjM^"*te39H- ?{F,P,i)P*&t)) 

= j&jfrrhl*%P*(&$ . (17) 



Aus der Gleichung (11) ist bekannt, dass in Gleichung (17) das Integral /dx P am (x, T) 
p*(x, T) = /dx F(x, T) P*(x, T) gleich dem europaischen Optionspreis P gu ist. Gilt nun 
P*(x, 0) = 5(x - S) als Anfangsbedingung des adjungierten Prozesses, erhalt man 
aus Gleichung (17) 

= pUS;0}+£ it jf'Mtm^^tmn^ (1 8) 

wodurch der fruhe Ausiibungsoptionspreis in Form eines Integrals uber den Fehl- 
preis-Quellterm Jausgedriickt wird, der durch den adjungierten Prozess P* gewichtet 

wird. 

Unter Verwendung der Gleichung (16), kann Gleichung (18) weiter vereinfacht wer- 
deh zu: 

jfc ^P^ + £m y| - Pfc4 (19) 

wobei n(t) das Gebiet der Ausubung zum Zeitpunkt t ist. 

Gleichung (16) ist hierbei ein Generalisierung der Integralgleichungen, die von Kim, 
Jacka, Carr, Jarrow und Myeni fur das Black-Scholes-Modell mit einem einzigen Ba- 
siswert hergeleitet worden sind. Diese Formel konnte beispielsweise - wie spater 
noch dargestellt wird - fur semi-analytische Methoden benutzt werden, urn eine Aus- 
iibungsgrenze zu finden. Es gilt hierbei, dass im Black-Scholes-Modell die Gleichung 
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(19) unter bestimmten Voraussetzungen die Zerlegungsformel ergibt, die in diesen 
Arbeiten hergeleitet wird. 

Die numerische Verwendung der APC-Methode in einer Finite-Differenzen-Methode 
wird im Folgenden erlautert. Es wird dabei von der Zerlegung in den Gleichungen 
(18) und (19) als ein Werkzeug ausgegangen, das die numerischen Methoden zur 
Losung des Problems der Schatzung von amerikanischen Optionen verbessert 

Es sind drei Eigenschaften, die diese Gleichung fur die numerische Bestimmung von 
Optionspreisen interessant machen: 

1. Der Quellterm ist nur in dem Gebiet der Ausubung von Null verschieden. In die- 
sem Gebiet ist auch die genaue Losung bekannt, da sie gleich der Payoff-Funktion F 
ist. Daher kann der Quellterm mit einer hohen Genauigkeit entweder numerisch oder 
analytisch berechnet werden. 

2. Wenn der adjungierte Prozess analytisch bestimmt werden kann, kann der Beitrag 
dieser Quellen zur Optionspreiskorrektur mit einer genauen Transferfunktion P*(S,T; 
S 0 ) zu S 0 propagiert werden. 

3. Bei der Monte-Carlo-Methode glattet das Doppelintegral Fluktuationen in den loka- 
len Preisschatzungen. 

Die Gleichungen (18) und (19) konnten zudem dazu verwendet werden, die Preis- 
schatzung einer einfachen expliziten Finite-Differenz-Methode zu verbessern. Zur 
Verdeutlichung wird ein Black-Scholes-Operator der Gleichung (8) verwendet. Das 
Einsetzen von 



- 15- 



mit a = - Vz (k - 1), p = - 1 A (k + 1f,k = 2r/<? transformiert die Black-Scholes- 
Gleichung in eine selbst-adjungierte Form 

fyp _ cft p 

welche numerisch sehr viel einfacher zu losen ist. Die Payoff-Fun ktion F(S) = max(E 
- S, 0) wird zu folgender Gleichung transformiert: 



. J 



ip,¥j = m^(eW-m _ ****** 0) (22) 

und die AnfangsbedingUngen sind 

p(x, 6) = f {x, 0) ^mh^®*^ 0 ~ fy. (23 ) 

Die Randbedingungen konnen nunmehr wie folgt hergeleitet werden: 

m = {24) 

Zur Verdeutlichung der APC-Methode im Finite-Differenz-Schema wird eine explizite 
Finite-Differenz-Diskretisierung der Gleichung (21) mit einem normalen Gitter { x 0 , x u 
x 2 , x N }> x n = x 0 + nSx und mit einem konstanten Zeitschritt St verwendet. Mit der 
Definition p™ :=p(x n ,mSt) ergibt sich die Diskretisierung fur eine europaische Put- 
Option aus 

Vernachlassigt man die Terme der Order O(St) und 0((Stf) erhalt man die Differenz- 
gleichung 

K +1 = *1%H + (> " »fe#***ft<i (27) 
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wobei 

-mm 



w (28) 



gilt. Es ist hierbei bekannt, dass diese Update-Methode stabil fur 0 < a < V2 ist. 

Fur eine amerikanische Option muss die Gleichung (27) durch folgende Gleichung 

ersetzt werden: 

maxfeti * (1 - + <mUrM^t (29) 

wobei f n m =f(x 0 +nSx 9 mSt) die diskretisierte Payoff-Funktion ist. Die Anfangs- und 
Grenzgleichungen sind durch die Gleichungen (23) bis (25) bestimnit. 

Aufgrund adjungierte Prozesskorrekturen im Finite-Differenz-Schema wird die Glei- 
chung (16) zur Berechnung des friihen Ausubungspreises als Zusatz zu der Finiten- 
Differenz-Methode verwendet. Gleichung (21) liegt in selbst-adjungierter Form vor, 
so dass 

•0- = 4 = 



m (30) 

gilt. Die Losung des in Gleichung (10) definierten adjungierten Problems ist 

fi?h r;*oM e*^N' (31) 

wobei x 0 := ln(S</£) und S 0 der Preis des zu ermittelnden Basiswerts ist, das heiBt in 
den meisten Fallen der derzeitige Preis. 

Der Ausubungspreis P am - P eu konnte dann durch die Aufsummierung von max(F- 
P am - LP am (dt, 0)) an jedem Gitterpunkt und bei jedem Zeitschritt berechnet werden. 

Fig. 1 zeigt die Abhangigkeit einer amerikanischen Option bei einem risikolosen 
Zinssatz r= 0,1, der Volatilitat cr= 0,4 und dem Ausubungspreis E= $ 10, wobei die 
Option in sechs Monaten ablauft und sie auf einem groben Gitter mit der expliziten 
Finite-Differenz-Methode (50 Gitterpunkte) berechnet worden ist. Fig. 1 zeigt ferner 
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auf dem gleichen Gitter die Abhangigkeit einer amerikanischen Option durch den 
adjungierten Prozesskorrekturterm, der zum europaischen Preis addiert wird. (50 
Gitterpunkte, APC-Methode), Fig. 1 zeigt zudem eine Referenzlosung, die durch ein 
Crank-Nicolson-Schema auf einem feinen Gitter berechnet ist, und den entsprechen- 
den Preis der europaischen Option. 

Die Ergebnisse fur eine Simulation auf einem groben Gitter sind in Fig. 1 gezeigt. 
Offensichtlich enthalt man erheblich bessere Schatzungen des amerikanischen Opti- 
onspreises, wenn der Preis durch das Addieren des fruhen Ausubungspreises, der 
mittels der Gleichung (18) berechnet wird, zu dem europaischen Optionspreis erfolgt, 
als bei der direkten Berechnung mittels der Finiten-Differenz-Methode, obwohl bei 
beiden Methoden das selben Gitter verwendet wird und die Summierung der gewich- 
teten Quellenterme nur einen geringen Rechenaufwand bedarf. Weiterhin ist ersicht- 
lich, dass die verknupfte Prozesskorrektur eine glatte Schatzung fiber dem gesamten 
Intervall ergibt, wahrend die Finite-Differenz-Methode zwischen den Gitterpunkten 
eine Interpolationsmethode verwendet. 

Urn die Verwendung der Gleichung (15) fur Monte-Carlo-Methoden - das heiBt die 
Schatzung von amerikanischen Optionen durch die stochastische Simulation - zu 
zeigen, wird sie auf die stochastische Simulation einer Mastergleichung angewendet. 
Dies ist ein numerischer Ansatz, der viele Anwendungen in der Physik hat, urn eine 
Mastergleichung zu losen, die eine Zeitentwicklung einer Wahrscheinlichkeit P n be- 
schreibt das System im Zustand n zu finden: 



wobei R nm die Ubergangsrate fur den Ubergang des Systems vom Zustand m in den 
Zustand n ist. Die Gleichung (32) muss hierbei als ein System von ordentlichen Diffe- 
rentialgleichungen fur alle Zustande gelesen werden. 

Wiederum wird die Anwendung der adjungierten Prozesskorrektur in einem Black- 
Scholes-Modell gezeigt. Auf jeden Fall ist dieser Ansatz nicht auf dieses Modell limi- 
tiert und Effekte wie die stochastische Volatilitat und Sprungprozesse konnen auch in 
einfacher Art und Weise in die Mastergleichung integriert werden. 




(32) 
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Trotz ihrer Flexibility und Einfachheit ist die stochastische Simulation keine Stan- 
dardmethode, die in der Finanzmathematik verwendet wird. Faller und Petruccione 
haben einen Mastergleichungsansatz fur die Optionspreisschatzung prasentiert, der 
sich allerdings in wichtigen Punkten von der hier gezeigten Methode unterscheidet, 
beispielsweise da diese konstante Zeitschritte verwenden. Der stochastische Simula- 
tionsansatz von Faller und Petruccione ist zwar schneller als die Standard-Monte- 
Carlo-Methoden, aber nicht so effektiv wie der hier gezeigte Ansatz. 

Nunmehr wird die Verwendung der Minimal-Prozess-Methode erlautert. Die Idee hin- 
^ \, ter der Minimal-Prozess-Methode ist es, einen einfachen Markov-Prozess zu erstel- 
len, dessen Zustandswahrscheinlichkeiten der Mastergleichung (32) gehorchen. 
Wenn eine Darstellung der Zustande gefunden ist, dann arbeitet der Algorithmus wie 
folgt: 

1. Auswahl eines Anfangszustands n gemaB den Anfangsbedingungen. 

2. Aufaddieren der Ubergangsraten fur alle Zustande m, die vom Zustand n in einem 
einzigen Ubergang erreicht werden konnen: R tot =^ m R mn - 

3. Vorschreiten in der Zeit urn ein Zeitinkrement x, das aus einer Exponentialvertei- 
lung erhalten wird und das Mittel r tot aufweist, das heiBt mit der Wahrscheinlichkeits- 
dichte tt(t)=R; 0 ) exp(-R tot r). 

4. Auswahl eines neuen Zustands m aus alien moglichen Zustanden mit der Wahr- 
scheinlichkeit R mr /R tot . 

5. Durchfuhren des Ubergangs zu dem ausgewahlten Zustand m, das heiBt es wird n 
= m gesetzt 

6. Gehe zu Schritt 2 bis die zunehmende Zeit groBer oder gleich der gewunschten 
Simulationszeit ist. 
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Der zeitaufwendigste Teil dieses Algorithmus ist der Schritt 4, die Auswahl des wirkli- 
c«hen Ubergangs, da die Nummer der Endzustande, die von m aus erreicht werden 
kann, sehr groB werden kann. Bei dem folgenden Beispiel ist es O(N), wobei N der 
Anzahl der verwendeten Gitterpunkte entspricht. Ein linearer Auswahlalgorithmus 
wurde demnach auch O(N) sein. Hier wird, um den Auswahlprozess zu beschleuni- 
gen, der ,,Dual-Logarithmic-Classes"-Ansatz von Fricke und Schnakenberg benutzt, 
bei dem alle moglichen Ubergange in Reaktionskategorien organisiert werden. In- 
nerhalb einer Kategorie variiert die Ubergangsrate hierbei nur hochstens um einen 
Faktor Zwei, so dass Ubergange sehr effizient durch ein von-Neumann- 
Verwerfungsschema ausgewahlt werden, wenn die Reaktionskategorie einmal ge- 
wahlt ist. Die Anzahl der Reaktionskategorien ist hierbei klein (die Anzahl von Okta- 
ven, die durch die minimale und maximale Ubergangsrate aufgespannt wird), so 
dass nur wenige Schritte benotigt werden, um eine Reaktionskategorie mit der 
Wahrscheinlichkeit R/R^ durch lineare Selektion auszuwahlen, wobei R c die Summe 
aller Ubergangsraten in der Kategorie c ist. 

Der Einfachheit halber konnte - als Diskretisierung der Black-Scholes-Gleichung 
durch eine Mastergleichung - im Folgenden wieder die Diffusionsgleichung verwen- 
det werden, das heiBt die in eine selbst-adjungierte Form transformierte Black- 
Scholes-Gleichung. 

Der Prozess wird unter folgenden Annahmen diskretisiert: 

1. Es wird eine aquidistante Diskretisierung des Logarithmus' des Preises des Ba- 
siswerts vorgenommen von x 0 = x min bis x N =x max m\t x t -x h1 = Ax fur /= 1,...N. 

2. Die Zustande sind fur jeden Raum [Xj, x M ) konstante Sprungfunktionen - wie in 
Fig. 2 dargestellt, welche die Diskretisierung des Preises des Basiswerts und die zu- 
gelasseneri Ubergange fur einen Zustand zeigt. 

3. Fur jeden Raum wird ein Sprung eines Quasi-Teilchens zu seinem nachsten linken 
oder rechten Nachbar zugelassen. Das Teilchen tragt einen test definierten Anteil Ap 
des Optionspreises, so dass beispielsweise nach einem Sprung von / zu / - 7, die 
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aktualisierten Preise p(i, t) = p(i, t - r) - Ap und p(i - 1, t) = p(i - 1, t- r) + Ap fur 
eine europaische Option und p(i, t) = max(p(i, t- r) - Ap, f(i, t)) und max(p(7- 1, t) = 
p(i- 1, t-r) + Ap, f(i ,t)) fur eine amerikanische Option betragen. 

Offensichtlich beschreibt dies einen zufalligen eindimensionalen Weg fur jedes Qua- 
si-Teilchen. Demnach unterliegt die Wahrscheinlichkeitsdichte einer Gruppe von 
Quasi-Teilchen einer Diffusionsgleichung mit der Grenze Ax ^ 0. In der selbst- 
adjungierten Form der Gleichung (21) ist der Diffusionskoeffizient gleich Eins, so 
dass die Ubergangsrate fur einen Sprungprozess von einem Quasi-Teilchen von /zu 
i± 1 unabhangig von /auf 

r = (33) 

gesetzt wird. 

Die Effizienz der stochastischen Simulation konnte ferner durch die Einfuhrung von 
Nettoubergangsraten verbessert werden. Die Anzahl von Quasi-Teilchen im Raum i 
zum Zeitpunkt fist N t = p(i, t) / Ap. Demnach ergibt sich aus der Mittlung aller Ober- 
gange, die am derzeitigen Zustand beginnen, dass der Nettotransfer der Quasi- 
Teilchen vom Knoten /zu i± 1 sich aus 

h " t o : msm*. (34) 

ergibt, wobeiAoo die durchschnittliche Gr63e des Zeitschritts ist und dies dadurch 
berucksichtigt wird, dass die Nettoubergangsrate fur jeden Sprung auf jeden Raum 
eingefuhrt wird. 

r^-t ij - | o ? Nr£W±i (35) 

Diese Nettoubergangsraten zu verwenden, hat zwei besondere Vorteile. Zum einen 
werden die Bewegungen der Belegungszahlen in den Raumen, das heiBt in den 
Preisschatzungen, stark reduziert. Zum anderen ist die Summe uber die Nettouber- 
gangsraten weitaus geringer als die Summe uber die Ubergangsraten der Quasi- 
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Teilchen. Da die typischen Zeitschritte umgekehrt proportional zu der Gesamtuber- 
gangsrate ist, bedeutet dies, dass die Simulation deutlich schneller beendet wird. 

Fig. 3 Der Preis einer amerikanischen Put-Option wird auf einen groben Gitter mittels 
einer kleinen Anzahl von Ubergangen (MC30) berechnet und durch Hinzufugen der 
adjungierten Prozesskorrektur in der gleichen Simulation der europaische Options- 
preis (APC) geschatzt. Die Volatilitat in diesem Beispiel ist a = 0,4, der Zinssatz r = 
0,1 und der Ausubungspreis E= $10. 

Als Beispiel fur ein numerisches Experiment zeigt Fig. 3 die Ergebnisse einer Schat- 
zung fur eine amerikanische Option, bei der ein sehr grobes Gitter mit 30 Gitterpunk- 
ten verwendet wird und 600 Ubergange in einer einzelnen Realisierung des sto- 
chastischen Prozesses durchgefuhrt werden. Die Optionspreise, die durch eine 
Ruckwartssimulation geschatzt werden, sind praktisch wertlos: bei manchen Gitter- 
punkten liegt der vorhergesagte amerikanische Preis unter dem europaischen Preis, 
was unmdglich ist. Die APC Methode basiert auf der gleichen Simulation und liefert 
dagegen bei nur gering erhohtem Rechenaufwand sehr gute Schatzungen des wirkli- 
chen amerikanischen Optionspreises. 

AbschlieBend und zusammenfassend ist festzustellen, dass der fruhe Ausubungs- 
preis fur eine amerikanische Option durch eine kompakte Formel dargestellt werden 
kahn, in dem von einem adjungierten Prozess Gebrauch gemacht wird, der sich in 
der ruckwartsgerichteten Zeitrichtung des Originalprozesses entwickelt, wobei die 
Entwicklung des entsprechenden europaischen Optionspreises bestimmt wird. Diese 
Formel generalisiert in erfinderischer Weise eine Integralgleichung, die irh Black- 
Scholes-Modell gultig ist und die in der Literatur dazuverwendet wird, semi- 
analytische Methoden zur Schatzung von amerikanischen Optionen abzuleiten. 

Es wird gezeigt, dass die Formel als Zusatz fur numerische Standardmethoden ver- 
wendet werden kann, urn entscheidend bessere Schatzungen des frCihen Aus- 
ubungspreises zu erhalten. Dies wird am Beispiel des Black-Scholes-Modells fur eine 
Finite-Differenz-Methode gezeigt. Um den Gebrauch der Formel fur Monte-Carlo- 
Methoden zu zeigen, wird eine stochastische Simulation, die auf einer Masterglei- 
chungsdarstellung des Prozesses basiert, eingefuhrt. Auch fur sehr grobe Annahe- 
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rungen mit schlechten Statistiken kann so eine hinreichend gute Schatzung des fru- 
hen Ausubungspreises erhalten werden. Der Grund fur die gute Leistung ist, dass es 
die Formel erlaubt, eine analytische oder semi-analytische Propagation des Quell- 
terms durchzufuhren, der in dem Ausubungsgebiet generiert wird. Demnach liegt die 
gr63te Ursache fur numerische Fehler in der Unsicherheit des Ortes der Ausubungs- 
grenze. 

Offensichtlich kann der adjungierte Prozess auch dazu benutzt werden, um neue 
semi-analytische und numerische Methoden zu entwickeln. Wenn man beispielswei- 
se eine Anfangsschatzung der Ausubungsgrenze hat, kann die stochastische Simu- 
lationsmethode derart verandert werden, dass ohne Diskretisierung in dem Aus- 
ubungsgebiet gearbeitet werden kann, da die Preise und damit die Ubergangsraten 
dort analytisch bekannt sind. In den folgenden Gebieten, konnen die Werte durch die 
Aufsummierung der adjungierten Prozesskorrekturen uber das Ausubungsgebiet be- 
rechnet werden, so dass bei einer selbstkonsistenten Berechnung, eine sehr grobe 
Diskretisierung um die Grenze ausreicht. In Kombination mit der Tatsache, dass 
durch einen adjungierten Prozess erhaltene Schatzungen auch bei sehr groben Git- 
tern verwendbar sind, folgt hieraus, dass die Methode auch sehr gut bei der Schat- 
zung von Optionen mit mehreren Basiswerten verwendet werden kann. 
SchlieRlich eroffnet der adjungierte Prozess eine besonders gute Moglichkeit, Preise, 
die auch groben Gittern berechnet sind, zu interpolieren and zu extrapolieren. 

In der Literatur ist folgende explizite Zerlegung des Preises einer amerikanischen 
Put-Option in den europaischen Preis und den fruhen Ausubungspreis im Black- 
Scholes-Modell bekannt: 

^imM = mm *m j*<r«m ^ MAl^zM^ i m (36) 

wobei B t der Ort der Ausubungsgrenze zum Zeitpunkt t ist und es gilt 

Dies ist ein spezieller Fall der Gleichung (19), denn bei L = L BS aus Gleichung (8) 
ergibt sich der adjungierte Operator aus 
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wie dies durch die partielle Integration in der bestimmenden Gleichung des adjun- 
gierten Operators gezeigt werden kann, siehe Gleichung (3). Der adjungierte Opera- 
tor ist hierbei nicht mit dem Black-Scholes-Operator identisch, hat aber einen ahnli- 
chen Aufbau. Das adjungierte Problem 

= ( 39 ) 

at 

kann analytisch durch die gleichen Methoden gelost werden, die zur Herleitung der 
Black-Scholes-Formel benutzt werden, wobei wesentlich einfachere Anfangsbedin- 
gungen verwendet werden konnen. Das Ergebnis ist: 



>i/%vmH ' (40) 
Die Payoff-Funktion einer amerikanischen Put-Option ist 

FW,$ = F0^= rnax|£?-4Q) (41) 



so dass innerhalb des Ausubungsgebiets gilt 



(42) 



Das Ausubungsgebiet ist in diesem Fall das Intervall Q(f) = [0, £J und aus Gleichung 
(19) folgt, dass 

f* * f Bi ^ rE , fef iMz fc ~ j£ggftt! \ (43) 

io /o % 2CT 2 t y/W&H 
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Durch die Substitution In S/S 0 und quadratische Erganzung im resultierenden Ex- 
ponentialargument kann dies in die Form der Gleichung (36) gebracht werden. 

Hinsichtlich weiterer vorteilhafter Ausgestaltungen und Weiterbildungen des erfin- 
dungsgemaRen Verfahrens wird zur Vermeidung von Wiederholungen auf die allge- 
meine Beschreibung sowie auf die beigefugten Patentanspruche verwiesen. 

SchlieBlich sei ausdrucklich darauf hingewiesen, dass die voranstehend beschrie- 
benen Ausfuhrungsbeispiele lediglich zur Erorterung der beanspruchten Lehre die- 
nen, diese jedoch nicht auf die Ausfuhrungsbeispiele einschranken. 



-25- 



Patentanspruche 

1. Verfahren zur Bestimmung der, insbesondere zukunftigen, Position eines Ob- 
jekts, insbesondere eines Teilchehs Oder dergleichen, wobei der Positionsverlauf 
eines Referenzobjekts bestimmt wird, wobei das Referenzobjekt einen Prozess 
durchlauft und wobei die Differenzen der Positionen des Objekts und des Differenz- 
objekts gebildet werden, 

dadurch gekennzeichnet, dass aus den Differenzen Fehlpositionsterme 
gebildet werden und dass die Fehlpositionsterme mit mindestens einem Korrektur- 
faktor gewichtet werden. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die mit dem Kor- 
rekturfaktor gewichteten Fehlpositionsterme aufaddiert werden. 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet der Korrekturfaktor 
durch eirien adjungierten Prozess berechnet wird. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass der adjungierte 
Prozess in entgegengesetzter Zeitrichtung des Prozesses des Referenzobjekts ver- 
lauft. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass zur Bestimmung 
der Position des Objekts die gewichteten und aufsummierten Fehlpositionsterme zu 
der Position des Referenzobjekts addiert wird. 
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Zusammenfassung 



Ein Verfahren zur Bestimmung der, insbesondere zukunftigen, Position eines Ob- 
jekts, insbesondere eines Teilchens Oder dergleichen, wobei der Positionsverlauf 
eines Referenzobjekts bestimmt wird, wobei das Referenzobjekt einen Prozess 
durchlauft und wobei die Differenzen der Positionen des Objekts und des Differenz- 
objekts gebildet werden ist im Hinblick auf eine gute Bestimmung der Position des 
Objekts bei gleichzeitig moglichst geringem technischem Aufwand derart ausgebil- 
det, dass aus den Differenzen Fehlpositionsterme gebildet werden und dass die 
Fehlpositionsterme mit mindestens einem Korrekturfaktor gewichtet werden. 



(Fig. 1) 
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